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A  computer  program,  labeled  SUMARY,  has  been  written  to  summarize  sound 
speed  values  computed  from  oceanographic  station  data.  The  program  determines 
the  following  values  at  each  of  40  standard  depths  from  selected  data  popula¬ 
tions:' 

mean  , 
variance 

standard  deviation 
minimum- 

10th  through  90th  percentiles, 
first  quantile' 
third  quantile 
max i mum  , 

number  of  observations  , 

Summary  data  is  available  on  punched  cards  for  146  of  240  ocean  data  populations 
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FOREWORD 


This  report  Is  one  of  a  series  describing  work  to  develop  summary  sound 
speed  profile  data  for  application  In  a  variety  of  underwater  acoustic  problems 
The  work  Is  under  joint  sponsorship  of  the  Long  Range  Acoustic  Propagation 
Project  (LRAPP)  and  the  Acoustic  Environment  Support  Detachment  (AESD)  of  the 
Office  of  Naval  Research  (ONR).  The  Fleet  Numerical  Weather  Central  (FNWC), 
Monterey,  California  is  assisting  In  this  project  by  furnishing  computer  time 
and  physical  support.  The  Naval  Undersea  Center  (NUC),  San  Diego,  California 
is  furnishing  a  research  oceanographer  to  perform  and  direct  the  work  at 
FNWC,  Monterey. 


INTRODUCTION 


This  report  Is  documentation  for  a  FORTRAN  Extended  computer  program  which 
performs  statistical  manipulations  and  computations  on  oceanographic  station 
data  organised  by  time-space  populations  contained  on  magnetic  tape.  The  program 
was  written  for  operation  on  the  CDC  6500  computer  system  at  the  Fleet  Numerical 
Weather  Central  (FNWC),  Monterey,  California, 

The  acronym  given  this  program  is  SUMARY.  The  purpose  for  SUMARY  is  to 
produce  [’•''intouts  of  ordered  sound  speed  data  at  40  standard  ocean  depths; 
select  eleven  statistical  percentiles  and  maximum  and  minimum  values  from  the 
ordered  set:;  compute  mean,  variance  and  standard  deviation  from  sound  speed 
data  at  each  of  the  40  depths;  and,  finally,  to  punch  the  selected  percentiles, 
maximum,  minimum,  mean,  variance,  standard  deviation  and  number  of  observations 
into  cards. 

The  40  standard  depths  used  for  SUMARY  are  shown  in  Table  I. 

The  current  report  is  the  third  in  a  series  describing  work  to  develop 
suaiary  sound  speed  data  by  specific  ocean  regions  for  application  in  a  variety 
of  underwater  sound  propage 'ion  problems.  References  1  and  2  are  other  reports 
in  this  series. 


INPUT 


Input  data  to  SUMARY  is  oceanographic  station  data  a:sembled  on  magnetic 
tape  by  a  program  labeled  TPC0NV.  These  data  are  sound  speed  (C),  temperature  (T), 
salinity  (S),  and  slgma-t  (o^)  at  40  standard  lev»c's.  The  data  are  grouped  by 
specific  ocean  regions  and  by  three  month  time  periods.  Appendix  A  is  a  layout 
of  the  input  data  tape  format.  NUC  Tech-ical  Note  1223  presents  a  description 
of  the  generating  program  TPC0NV. 


OUTPUT  DATA 


SUMARY  produces  a  series  of  printout  listings  and  a  deck  of  data  cards. 

The  computer  prints  out  40  tables  of  ordered  sound  speed  values.  One  for 
each  of  the  40  standard  depths  shown  in  Table  I.  Tne  maximum  value  heads  each 
list,  the  minimum  value  is  at  the  end.  The  tables  contain  twenty  decreasing  values 
per  line.  The  number  of  data  values  determines  the  number  of  iines.  Table  II 
is  a  typical  ordered  data  list  for  scund  speed  values  at  depth  1 ,  or  0  meters. 

The  printed  values  are  sound  speed  in  meters  per  second  minus  1000. 
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TABLE  I 

STANDARD  DEPTHS  AND  ASSOCIATED  DEPTH  NUMBERS  USED  FOR  SUMARY 


*TH  NO. 

DEPTH 

1 

0 

2 

10 

3 

20 

4 

25 

5 

30 

6 

35 

7 

50. 

8 

75 

9 

100 

10 

125 

11 

150 

12 

200 

13 

250 

14 

300 

15 

400 

16 

500 

17 

600 

18 

700 

19 

800 

20 

900 

21 

1000 

22 

1100 

23 

1200 

24 

1250 

25 

1300 

26 

1400 

27 

1500 

28 

1750 

29 

2000 

30 

2500 

31 

3000 

32 

4000 

33 

5C00 

34 

6000 

35 

7000 

36 

8000 

37 

9000 

38 

10000 

39 

11000 

40 

12000 
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A  second  type  of  output  Is  the  deck  of  data  cards.  This  outout  card  deck 
contains  80  cards,  two  cards  of  Information  for  each  of  the  40  standard  depths. 

The  first  card  of  each  pair  contains  maximum,  10,  20,  30,  40,  50,  50,  70,  80, 
and  90th  percentiles,  minimum,  first  quantile,  and  third  quantile  sound  speed 
values  at  a  specific  depth.  This  card  also  contains  the  depth  number  (see  Table  I), 
the  number  of  sound  speed  observations  In  the  population  and  the  specific  label 
given  the  data  group.  The  second  card  for  a  specific  depth  contains  standard 
deviation,  variance,  number  of  observations  in  decimal  form,  mean  value,  depth 
number,  number  of  observations  in  Integer  form  and  the  date  group  label.  The 
program  inserts  a  flag  card  containing  80  columns  of  nine  (9/  to  separate  data 
card  groups.  Appendix  B  Is  a  description  of  the  output  data  card  format. 

Appendix  C  is  a  typical  listing  of  an  output  data  card  deck. 
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GENERAL  PROGRAM  OPERATION 


SUMARY  was  written  to  operate  sper < f leal ly  on  sound  speed  data.  A  simple 
program  modification  will  permit  operation  on  any  one  of  the  other,  T,  S,  or 
ot  data  populations. 

SUMARY  will  process  one  or  more  populations  of  sound  speed  data  in  a 
single  run.  A  series  of  data  cards  directs  the  computer  in  selecting  the  appro¬ 
priate  data  groups,  or  populations,  from  the  input  tape  for  processing.  Appendix 
D  is  a  description  of  the  data  cards  format. 

The  computer  central  memory  Is  insufficient  to  contain  an  entire  potential 
data  population  at  each  of  the  40  standard  depths.  For  this  reason  the  program 
reads  data  from  the  input  tape  and  fills  a  holding  buffer  in  central  memory 
then  outputs  this  buffer  to  disk  storage.  The  program  reads  blocks  of  data, 
lords  the  holding  buffer  up  to  500  observations  at  each  of  the  40  standard  depths, 
th»n  transfers  the  full  buffer  to  disk  storage  before  accepting  more  data.  This 
process  continues  until  a  complete  data  population  has  been  transferred  from 
tape  to  disk.  Upon  completing  the  Input  transfer  of  data  the  program  begins  to 
operate  on  the  data  stored  on  disk.  All  data  for  one  depth  Is  brought  from  disk 
back  into  central  memory.  Once  the  entire  data  population  from  one  depth  is  in 
central  memory  the  program  orders  the  values,  computes  mean,  variance  and  stan- 
datd  deviation  and  selects  the  minimum,  maximum  and  eleven  percentile  values. 

The  ordered  list  Is  then  printed  out  and  the  appropriate  data  cards  punched. 

This  procedure  continues  through  the  40  standard  depths.  The  program  then 
determines  if  another  data  group  Is  to  be  processed  or  If  the  job  is  to  end. 

This  process  continues  until  a  data  card  directs  the  computer  to  end  the  job. 


DATA  CARDS 


The  series  of  cards  which  control  the  selection  of  data  groups  for  pro¬ 
cessing  include  one  or  more  data  cards  and  a  terminating  flag  card. 

Each  data  card  contains  three  fields.  The  first  field  is  found  in  col¬ 
umns  1  through  5.  The  computer  will  pass  over  the  number  of  data  groups  on 
the  input  tape  that  is  punched  into  this  field.  Processing  will  begin  on  the 
data  group  number  following  the  figure  in  this  first  field.  The  second  field 
is  found  In  card  columns  6  through  10.  This  figure  instructs  the  computer  in 
the  number  of  data  groups  to  process  under  the  current  data  card  control.  The 
computer  will  process  the  number  of  consecutive  data  groups  indicated  by  the 
number  in  this  second  field.  The  last  field  is  found  In  columns  11  through  15. 
This  is  an  optional  field.  This  field  contains  the  number  of  E0C  cn  the  input 
data  tape.  If  field  three  Is  blank  the  number  of  EOF  is  assumed  to  be  1.  A 
separate  data  card  must  appear  in  the  data  card  series  for  each  non-consecutive 
data  group  to  be  processed. 

A  flag  card  containing  nines  (9* s )  in  columns  1  through  80  signals  the 
computer  to  end  the  job. 
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PROGRAM  FORM 


The  basic  program,  SUMARY,  Is  In  the  form  of  a  400  card  source  deck.  The 
general  program  operation  Is  run  from  this  source  deck. 


PROGRAM  OPERATING  PROCEDURE 

The  procedure  for  obtaining  a  production  run  of  SUMARY  Is  to  follow  this 
sample  procedure. 

1.  Select  the  Input  data  tape. 

2.  Punch  the  appropriate  system  control  card  to  request  the 
selected  data  tape. 

3.  Determine  the  data  groups  to  be  processed  from  the  Input  tape. 

4.  Punch  the  appropriate  data  and  flag  cards. 

5.  Insert  the  data  cards  at  the  end  of  the  SUMARY  source  deck. 

6.  Submit  for  production  run. 


CONCLUSION 


SUMARY  Is  a  FORTRAN  Extended  computer  program  designed  to  (1)  generate 
ordered  lists  of  sound  speed  data  at  each  of  40  standard  levels;  (2)  select 
maximum,  minimum  and  eleven  percentiles  from  each  ordered  list;  (3)  compute 
mean,  variance,  standard  deviation  and  number  of  observations  at  each  depth; 
and  (4)  punch  the  selected  and  computed  values  Into  a  deck  of  data  cards. 

The  program  was  written  to  automatically  select  and  process  sound  speed 
data.  A  simple  program  modification  will  permit  the  selection  and  proceslsng 
of  temperature,  salinity  or  slgma-t  data. 

Production  runs  for  sound  speed  data  have  been  run  for  many  ocean  areas. 
Figure  1  Is  a  chart  showing  60  ocean  regions.  NUC  Technical  Note  1223, 
Reference  2,  discusses  the  origin  of  Figure  1.  Table  1111$  a  summary  of  the 
areas  shown  In  Figure  1  for  which  outputs  from  SUMARY  are  available. 

A  program  listing  Is  given  In  Appendix  E. 
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FIGURE  1.  Marsden  Square  Chart  Showing  Sixty  Designated  Ocean  Regions 
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AREA 

101 

102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 

Note: 

X  = 

G  = 


TABLE  HI 


SUMMARY  OF  DATA  PROCESSED  THROUGH  SUMARY 
SEASON  SEASON 

1  2  3  4  AREA  1  2 


X  X  X  X 

X  X  X  X 

X  C  X  X 

X  X  X  X 

X  X  X  X 

X  X  X  X 

X  X  X  X 

O  G  O  c 

X  G  X  X 

X  X  X  X 

X  X  X  X 

x  x  O  C 

G  X  G  G 


201  X  X 

202  X  X 

203  X  X 

204  X  X 

205  X  X 

206  X  X 

207  O  X 

208  X  X 

209  X  X 

210  X  X 

211  X  X 

212  X  X 

213  X  X 


X  X  X  X 
X  X  X  X 
X  X  X  X 
X  X  X  X 
X  X  X  X 
X  X  X  X 
X  X  O  X 


214  X  X 

215  X  X 

216  X  X 

217  X  X 

218  O  X 

219  X  X 

220  X  X 

221  C  C 


Data  processed  through  SUMARY 
Data  not  processed  through  SUMARY 
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APPENDIX  A 


REQUIRED  MAGNETIC  TAPE  FORMAT  FOR  INPUT  TO  SUMARY 


*am, 


TPCONV  OUTPUT  DATA  TAPE  FORMAT 


WORD 

DATA 

BITS 

BIT  POSITIONS 

1 

Label  for  the  area 

15 

(46-60) 

Blank  " 

5 

(41-45) 

Standard  depth  Indicators 

40 

(1-40) 

2 

Latl tude 

19 

(42-60) 

Longi tude 

20 

(22-41 ) 

Day 

10 

(12-21) 

Month 

4 

(8-11) 

Year 

7 

(1-7) 

3 

Marsden  Square 

15 

(46-60) 

Layer  depth 

15 

(31-45) 

Layer  sound  speed 

15 

(16-30) 

Blank 

9 

(7-15) 

No.  of  standard  levels 

6 

0-6) 

4 

In-layer  gradient 

15 

(46-60) 

Below-layer  gradient 

15 

(31-45) 

Axis  depth 

15 

(16-30) 

Axis  sound  speed 

15 

(1  15) 

5 

Sound  speed 

15 

(46-60) 

Temperature 

15 

(31-45) 

Salinity 

15 

(16-30) 

Sigma-t 

15 

0-15) 

Repeat  word  5  for  each  recorded  standard  level. 

Standard  recorded  depths  In  meters:  0,  10,  20,  25,  30,  35,  50,  75,  100,  125, 
150,  200,  250,  300,  400.  500,  600,  700,  800,  900,  1000,  1100,  1200,  1250,  1300, 
1400,  1500,  1750,  2000,  2500,  3000,  4000,  5000,  6000,  7000,  8000,  9000,  10000, 
11000,  12000. 

Start  each  new  station  with  a  new  word  1 . 


Maximum  of  511  words  per  record. 

Three  special  records  are  placed  at  the  end  of  each  data  group. 


FORMAT  FOR  THREE  SPECIAL  records  following  a  data  assembly  group 
Record  1  -  5  words 


WORD 

DATA 

BITS 

BIT  POSITIONS 

1 

Area  label 

15 

(46-60) 

Blank 

45 

(1-45) 

2 

No.  records  output 

30 

(1-30) 

Blank 

30 

(31-60) 

3 

No.  stations  accepted 

30 

(1-30) 

Blank 

30 

(31-60) 

4 

Length  of  output  data 

30 

(1-30) 

Blank 

30 

(31-60) 

5 

Zeros 

60 

(1-60) 

Record 

c  -  5000  words 

This  record  contains  a  geographical  distribution  map  of  data  assembled 
In  the  current  data  group.  The  record  Is  written  as  MAP  (Latitude,  Longi¬ 
tude),  where  latitude  goes  from  1  to  50  and  longitude  from  1  to  100. 

Record  3  -  1200  words 

This  record  contains  a  time  distribution  by  months  and  years  table  of 
data  assembled  In  the  current  data  group.  The  record  Is  written  as  MYCNT 
(Year,  Month),  where  year  goes  from  1  to  100  and  month  from  1  to  12. 

The  three  5,  500,  and  1200  special  word  records  separate  data  groups. 

One  EOF  occurs  at  the  end  of  data  on  the  tape. 
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APPENDIX  B 


OUTPUT  DATA  CARD  FORMAT 


SUMARY  OUTPUT  DATA  CARD  FORMAT 


CARD  NO. 
1 


3 

4 


CARD  COLUMNS 

DATA 

2-5 

Maximum 

7-10 

10th  percentile 

12-15 

20th  percentile 

17-20 

30th  percentile 

22-25 

40th  percentile 

27-30 

5Cth  percentile 

32-35 

60th  percentile 

37-40 

70th  percentile 

42-45 

80th  percentile 

47-50 

90th  percentile 

52-55 

Minimum 

57-60 

1st  quartile 

62-65 

3rd  quartile 

69-70 

Depth  number (s^e  Note  4) 

71-75 

Number  of  observations 

77-80 

Numerical  area  label 

1-10 

Standard  deviation  -  decimal 

11-20 

Variance  -  decimal 

21r30 

Number  of  observations  -  decimal 

31-40 

Mean  sound  speed  minus  * 
1000.0000  decimal 

55-56 

Depth  number  -  integer 

58-62 

Number  of  observations  -  integer 

66-68 

Numerical  area  label 

1-80 

Repeat  of  Card  1  format  for 
depth  2 

1-80 

Repeat  of  Card  2  format  for 
depth  2 

B-2 


■*  t 


79 

1-80 

Repeat  of  Card  1  format  for 

depth  40 

80 

1-80 

Repeat  of  Card  2  format  for 

depth  40 

81 

1-80 

99* *.9  mandatory  flag 

NOTES: 


1.  There  are  two  cards  for  each  of  40  depths. 

2.  Card  formats  1  and  2  are  repeated  for  each  of  the  40  depths. 

3.  Depth  number  Is  the  sequence  number  of  depth  value  in  the 
following  list. 

0,  10,  20,  25,  30,  35,  50,  75,  100,  125,  150,  175,  200,  250, 
300,  400,  500,  600,  700,  800,  900,  1000,  1100,  1200,  1250, 
1300,  1400,  1500,  1750,  2000,  2500,  3000,  4000,  5000,  6000, 
7000,  8000,  9000,  10000,  11000,  12000. 

4.  All  percentile  sound  speed  values  (C)  of  card  &  type  are 
recorded  (10C -10000). 

5.  All  depth  values  are  in  meters. 

6.  All  sound  speed  values  are  in  meters  per  second. 


APPENDIX  C 


TYPICAL  OUTPUT  CARD  DECK  LISTING 
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•* .  o6 16 
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O.f.OOO 
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*■'  8  117  210* 

8316  8260  8998  5342  9 

r  '  9  US  2104 
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10  118  2104 

8295  ?235  5357  5310  11 

’  11  118  2104 

8250  8202  83*7  5270  12 

12  118  2104 

8193  5101  5257  5213  13 

13  115  2164 

8120  5047  5206  5186  14 

14  US  2104 

4991  4931  5102  5091  15 

15  115  2104 

49o7'4072  4995  4937  16 

16  113  2104 

4066  4645  4923  4078  17 

17  113  2104 

4043  4631  4005  4051  10 

16  112  2104 
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19  106  2104 
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20  100  2164 

4025  402a  4041  4001  21 

*1  73  2184 
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22  00  2104 

403^  4831  40*9  4098  23 

23  17  2104 

9  0  0  0  .2* 

24  0  0)64 

’845  4043  4053  4000  gB 
1  **  .kv  25  ,  7  2104 

*&?1t  #rw*> 


wr&t 


0  49 

xfoO 

3  *fe.5 

92* 

m 

0^14996.499^,* 

r** 


9  *104 
117  2104 
U?  «|B* 
115  2104 
115  2164 
115  2104 
115  2104 
115  2104 
US  2184 
115  2184 
113  2184 
113  2184 
112  2104 
106  2184 
100  2104 
73  2104 
60  2184 
17  2104 
0  2184 
»7  2104 
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APPENDIX  D 


DATA  CARDS  FORMAT 


D-l 


t/siWwrtj  .^^MSaaS.  i?yJ.*MSSi&i&&miiSi 


SUMARY  DATA  CARD  FORMAT 


CARD  COLUMNS  LEGAL  CHARACTERS 


NOTES 


1-5 


Integers 


No.  of  data  groups  to 
skip  on  Input  data  tape 
before  reaching  the  first 
data  group  to  be  processed. 


6-10 


Integers  No.  of  consecutive  data 

groups  to  process 


11-15 


Integers  or 
blank 


No.  of  files  on  Input  data 
tape. 

Blank  Implies  1  file. 


Additional  data  cards,  one  for  each  non-adjacent  continuous  sequence  of  data 
groups. 


1-80 


999... 9 


Mandatory  flag  following 
last  data  card. 
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SUKARY  PROGRAM  LiSTING 
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131  SUMC=C. 
SCS?=u. 
XNO~0« 
VAR=0. 
SISMA=0 

crs=o. 

PINC=G. 

XN=P. 

f’Tl=0. 
IPTi=C. 
PT2=r. . 
PI  3=0. 
PCTL=0. 


FEkINl  2 

CALL  IIuf.H(CP»9) 

F  iT  1  =•_  . 

IFaTJ=U 

FGT2=C. 

FGT3=&. 

GTL1=0. 
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DEPARTMENT  OF  THE  NAVY 

OFFICE  OF  NAVAL  RESEARCH 
875  NORTH  RANDOLPH  STREET 
SUITE  1425 

ARLINGTON  VA  22203-1995 


IN  REPLY  REFER  TO: 


5510/1 

Ser  32 1OA/0 11/06 
31  Jan  06 


MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1)  List  of  DECLASSIFIED  LRAPP  Documents 

1 .  In  accordance  with  reference  (a),  a  declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


^ ^ 
BRIAN  LINK 
By  direction 


Subj :  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

DISTRIBUTION  LIST: 

NAVOCEANO  (Code  N121LC  -  Jaime  Ratliff) 

NRL  Washington  (Code  5596.3  -  Mary  Templeman) 

PEO  LMW  Det  San  Diego  (PMS  181) 

DTIC-OCQ  (Larry  Downing) 

ARL,  U  of  Texas 

Blue  Sea  Corporation  (Dr.Roy  Gaul) 

ONR  32B  (CAPT  Paul  Stewart) 

ONR  3210A  (Dr.  Ellen  Livingston) 

APL,  U  of  Washington 

APL,  Johns  Hopkins  University 

ARL,  Penn  State  University 

MPL  of  Scripps  Institution  of  Oceanography 

WHOI 

NAVSEA 

NAVAIR 

NUWC 

SAIC 


Declassified  LRAPP  Documents 


